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We investigate theoretically properties of two-dimensional topological insulator constrictions both 
in the integer and fractional regimes. In the presence of a perpedicular magnetic field, the constric¬ 
tion functions as a spin filter with near-perfect efficiency and can be switched by electric fields only. 
Domain walls between different topological phases can be created in the constriction as an interface 
between tunneling, magnetic fields, charge density wave, or electron-electron interactions dominated 
regions. These domain walls host non-Abelian bound states with fractional charge and spin and 
result in degenerate ground states with parafermions. If a proximity gap is induced bound states 
give rise to an exotic Josephson current with Syr-peridiodicity. 
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Introduction. The field of topological properties in con¬ 
densed matter systems has been rapidly growing over 
the past decade. In particular, the topics of topologi¬ 
cal insulators (TIs) [THH] and exotic bound states with 
non-Abelian statistics have attracted a lot of attention 
theoretically and experimentally [T5ll46j . Of special in¬ 
terest are also quantum effects arising from geometric 
confinement such as topological insulator constrictions 
(TICs) [T7H5^ , where the edge modes get coupled by tun¬ 
neling and a gap is opened in the energy spectrum, see 
Fig.lll In general, such a coupling is not desirable since 
typically it leads to a suppression of topological proper¬ 
ties mm- However, we will find that, quite surprisingly, 
in the presence of additional mode-mixing perturbations 
such as magnetic fields, superconductivity, and interac¬ 
tion effects, the tunneling does not necessarily destroy 
all such properties. Instead, different topological phases 
can emerge that give rise to exotic phenomena such as 
fractional fermions, parafermions, and exotic supercon¬ 
ductivity where Cooper pairs themselves get paired. 

First, we consider constrictions with extended edge 
modes and show that in the presence of tunneling the 
TIC can be tuned between insulating, propagating, and 
spin filtering regimes by electric fields only, making such 
TICs attractive candidates for spintronics applications. 

Second, we focus on localized modes. Here, we iden¬ 
tify competing mechanisms that generate gaps in the 
spectrum, arising from magnetic fields, tunneling be¬ 
tween edges, periodic modulations of the chemical poten¬ 
tial, proximity effects, and electron-electron interactions. 
We show that there are zero-energy bound states at the 
interfaces between two phases controlled by competing 
gap mechanisms. These states are fractional fermions 
of the Jackiw-Rebbi type. If fractional TIs with frac¬ 
tional charge e/m are considered, the ground state is 
m-fold degenerate and the resulting bound states are 
Zrn-parafermions. The superconductivity that could be 
induced by proximity effect at such constrictions corre¬ 
sponds to the coherent tunneling of two Cooper pairs and 
results in an unusual Syr-periodic Josephson current. 


TIC Model. We consider a constriction created in a 
two-dimensional TI, see Fig. Upper and lower edges 
of the TIC hosting helical states of opposite helicities 
are brought close to each other and, as a result, cou¬ 
ple via tunneling. The edges of the TIC are labeled 
by the index r, where r = I (r = —I) corresponds to 
the upper (lower) edge. The helical edge states of the 
TIC have a linear energy dispersion. The correspond¬ 
ing kinetic part of the Hamiltonian is given by Hj-in = 
—ihvp '^^{R\.dxRT — LXdxLr), where vp is the Fermi ve¬ 
locity. The operator Rt{x) [L,-(x)] is the annihilation op¬ 
erator acting on the right-propagating (left-propagating) 
electron located at point x of the TIC edge r. We note 
here that the two pairs of helical edge states possess oppo¬ 
site helicities, i.e., right-propagating (left-propagating) 
electrons at the upper (lower) edge are spin-up electrons 
and left-propagating (right-propagating) electrons at the 
upper (lower) edge are spin-down electrons, see Fig. 

We point out that the same setup could be assembled 
by bringing close to each other two TI samples [IS] or in 
the framework of strip of stripes models |S3| . The latter is 
especially important for the fractional regime [SSHSS] as it 
allows one to design fractional TIs for an array of coupled 



FIG. 1. A sketch of a constriction formed in a topological in¬ 
sulator. Spin up (down) edge states are shown in red (blue). 
The constriction could be either doped with magnetic impuri¬ 
ties or subjected to a magnetic field. The pairs of helical edge 
states are coupled by tunneling which results in the opening 
of an energy gap at zero momentum. A magnetic field applied 
in the plane of the TI couples edge states with opposite spins 
and also opens a gap in the spectrum. 
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one-dimensional channels with spin-orbit interaction |53j . 

The tunneling in the TIC is assumed to be spin con¬ 
serving and described by Htun = + H.c., 

where t is the tunneling matrix element between two 
edges of the TI and the operator is the electron 
annihilation operator at position x of the edge state r. 
In what follows, we use the fact that the fast oscillat¬ 
ing part of the wavefunction is given by where 

kp is the Fermi wavevector set by the chemical potential 
/i = hvpkp. Keeping only slowly varying terms in the 
Hamiltonian [211 [MM5] , we arrive at 

Htun = t{R\Li + l\Ri + H.C.). (1) 

The TI surface could be subjected to a magnetic field 
or doped with magnetic impurities producing a local 
effective magnetic field. The Hamiltonian is given by 
Hz = where the unit vec¬ 

tor n points along the field and cr is a vector composed 
of Pauli matrices acting on the electron spin. For fields 
along the spin quantization axis of the edge states which 
is chosen, say, in the ^-direction, the corresponding ef¬ 
fective Zeeman term is given in terms of right and left 
movers as 

H, = - l\Li - r\Ri + l\Li + H.C.), ( 2 ) 

where is the coupling constant either determined by 
the Zeeman energy or by the strength of exchange inter¬ 
action. If A^ is generated by a magnetic field with 
the vector potential A = B^yx applied perpendicular 
to the TI plane, then the wavevector k gets shifted to 
k — {e/hc)Ax, accounting for orbital effects of the mag¬ 
netic field m- If the upper (lower) edge state is at j/ = d 
{y = —d), the shift is given by —T{e/Hc)d. Interestingly, 
for TI edge states, the orbital and spin contributions add 
up to /S.Z — e-B^dvp/ c. However, for typical TIC sizes one 
can neglect the orbital part. 

The magnetic field applied perpendicular to the spin 
polarization axis, say, in the x direction, results in the 
Hamiltonian 

H, = A,(4Li + r\Li + H.C.), (3) 

where A^, is the strength of coupling in the x direction. 

Spin filter effect. In the presence of both tunneling 
and magnetic fields, the total Hamiltonian is given by 
H = Hkin + Htun + Hz -l- Hx and can be rewritten in the 
basis (i?i, Li, i?i, Li) in terms of Pauli matrices as 

■H = UFPP3 + A,rPl + A^P3T3 -I-tpiTi, (4) 

where p = —ihdx is the momentum operator and, for sim¬ 
plicity, we assume that A^,, A^, and t are non-negative 
if not specified otherwise. The Pauli operators pi act in 
right/left mover space. The energy spectrum is given by 

={hvpkf + + 

± 2v/A2t2 + [{fivpky + ( 5 ) 


We are interested in the regime Az > t > Ax- First, 
we notice that the two Dirac cones are shifted by the 
perpendicular magnetic field by kz = Az/hvp to the left 
(right) for the upper (lower) edge. The tunneling opens 
a gap at zero momentum fc = 0 of the size Ak=o = 2t, see 
Fig.@ while the magnetic field in the x direction opens 
a gap at finite momentum k = ±kz and at zero energy, 
given by Ak=±k, = 2Ax, see Fig. |2] 

The described setup can be used as a spin filter con¬ 
trolled purely by electric gates. For example, if A^, = 0, 
the spin projection on the z axis, Sz, is a good quan¬ 
tum number and all modes are spin polarized. If the 
chemical potential lies in the electron (hole) part of the 
spectrum, p e (A^ -t,Az+ 1) [p € {-Az - 1, -Az + 1)], 
only the spin down (spin up) component can propagate 
through the TIC, see Fig. However, due to the tun¬ 
neling there is leakage from the upper to the lower edge 
such that the probability to stay in the upper edge is 
given by {tz)\s^ ~ 1 — {t/^Azff. If p is tuned close to 
zero, \p\ < Ax, the system is fully insulating. For other 
values of p both spin components can propagate. 

If Ax 0, the propagating modes are no longer per¬ 
fectly spin-polarized, however, deviations are small in the 
ratio AxjAz- The chemical potential should be tuned 
into the window of /r S A'l + A^ — t, ^/A^ -|- A^ -|- t) 
[p G i-ff'A^ + Al - t, -^/Al+Al + t)] for the spin 
down (up) dominated propagation, see Fig. The efh- 
ciency of the spin filter is characterized by the probability 
to keep an initial spin polarization and to stay at the ini¬ 
tial edge, which is given by 

{tz) Is. « 1 - (t/AAzf - {Ax/2Azf. (6) 



FIG. 2. The energy spectrum of the TIC in the presence 
of magnetic fields (or exchange interactions). A magnetic 
field in ^-direction, the spin quantization axis of the TIC, 
shifts the Dirac cone corresponding to the upper (lower) edge, 
T = 1 (t = 1), to the left (right) such that the Dirac point is 
located now at wavevector —kz (kz). The tunneling between 
edges opens a gap of size 2t at fc = 0. A magnetic field in 
X- direction, i.e., perpendicular to the spin quantization axis, 
opens a gap of the size 2Ax at zero energy. If the chemical 
potential pi {pi) is tuned inside the gap opened by tunneling 
at fc = 0, only spin-polarized modes with spin up (spin down) 
propagate through the TIC. 
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By changing the position of the chemical potential, i.e., 
by applying electric fields, one can tune the TIC into dif¬ 
ferent spin filtering regimes. This provides a substantial 
advantage over spin filters tuned by magnetic fields which 
are difficult to switch fast and locally. 

Bound states at the tunneling-magnetic field interface. 
The TIC with spectrum Eq. ([^ not only allows one to re¬ 
alize spin filtering but also to trap bound states that are 
localized at the interface between tunneling- and mag¬ 
netic field-dominated regions. While the energy branch 
is always gapped, the branch is gapless at zero 
momentum if = A^-|-A^. At other values is gapped 
unless Ax = 0. We note that, generally, the interface sep¬ 
arating two regions that are characterized by opposite 
signs of the expression — A^ — A^ hosts zero-energy 
bound states. As an example, we consider an interface at 
the left end of the TIC {x = 0) specified by t = 0 for x < 0 
and by t > A^, with = 0 for x > 0. This interface 
hosts a zero-energy bound state with wavefunction of the 
form <i)(x) = (/, —if, f*,if*)'^ with /(x) = i9{x)e~^^^* -\- 
i0{—x)e''^‘^e~^/^^, where the localization lengths are de¬ 
fined as = hvF/{t— Aa;) and = hvp/Ax- An analo¬ 
gous bound state occurs also at the right end of the con¬ 
striction. These zero-energy bound states are examples 
of fractional fermions of the Jackiw-Rebbi type [67H72] 
and possess non-Abelian braiding statistics m- 

In passing we note that, alternatively, degenerate 
bound states even occur for t = A^ = 0, namely in the 
presence of a magnetic domain wall separating two do¬ 
mains with Aa;(x) = 9(x)Ax — 9{—x)Ax. Such an in¬ 
terface hosts a zero-energy bound state at each of the 
two TI edge states. The corresponding wavefunction is 
<i)(x) = (f, where = hvp/Ax. We note 

that as the gap closes twice (at the upper and at the 
lower edge), the twofold degeneracy is not protected and 
states split away from zero if the tunneling is included. If 
the magnetization rotates not exactly by tt but by some 
finite angle y, the bound state moves away from zero en¬ 
ergy, E = — Aa;C0s(x/2) for 0 < y < 27r. The domain 
wall localizes the charge e/2 only for y = tt, which brings 
us back to the fractional fermions of the Jackiw-Rebbi 
type 

Bound states at charge density wave - magnetic field 
interface. An alternative way to generate bound states 
is to allow for modulations of the chemical potential fj-mod 
with the period of 2kz, fimod = 2Sficos{2kzx-\-(j)), where 
26fj, is the amplitude of modulations and is the phase 
at X = 0. This creates a charge density wave (CDW) 
that opens a gap around the Fermi points. This setup 
works in the spin-filtering regime A^ 3> t,Ax. Indeed, 
the corresponding Hamiltonian is given by 

Hmod = + H.c.), (7) 

where the coupling amplitude is found in second order 
perturbation theory as fi ^ tS^/A^. The energy spec¬ 
trum then becomes + (A^, ±/l)^. The bulk 


gap closes if Ax = fi, indicating the topological phase 
transition. Thus, we can construct an interface between 
the magnetic field dominated region with /i = 0 (x < 0 ) 
and the CDW dominated region with A^, = 0 (x > 0). 
Again, such an interface hosts a zero-energy bound state 
with wavefunction of the form <i)(x) = {f,if,f*,—if*)^ 
with /(x) = -I- 0(—x)e®/^<] in 

the basis of ('kn,'ku,'kji, d>n)- Here, the localization 
lengths are defined as ^> = hvp/jl and ^< = hvp/Ax. 

Fractional bound states at the charge density wave - 
magnetic field interface. Next, we consider helical edges 
of a fractional TI constriction with elementary excita¬ 
tions of charge elm defined in terms of chiral bosonic 
fields (frn [73]. First, we note that also in this regime the 
system could be operated as a spin filter for quasiparti¬ 
cles. Second, we focus on properties of domain walls in 
this system. The electron operators are then rewritten as 
and To satisfy the anticom¬ 

mutation relations between original fermionic operators, 
we work with the following non-zero commutators for the 
bosonic fields, 

%1TT 

[4>rn{.x),4>r'n'{x')] = — Jnn'Sgn(x - x'). (8) 

All other commutators are assumed to vanish. To pro¬ 
ceed, we bosonize the magnetic field Hamiltonian ELx [see 
Eq. ([^] as Fix = 2Aa;( cos[m(0ii — -I- cos[to((/ii — 

/'ll)]) the CDW Hamiltonian Hmod [see Eq. Q] as 
Hmod = 2/l(cos[m((/ii-(/n-/’)]+cos[m((/ii-(/ii-k/')]). 
In a next step, we express the chiral fields firn in terms 
of their conjugated fip and 9p fields [p = ±1) defined as 
4>rn = {4)i/m-\-r9i/m-\-n{(l>i-\-r9i)\/2. The commutation 
relations between the newly introduced fields are given 
by [^i(x), 0 i(x')] = inm sgn(x —x') and [(f)i{x),9i{x')\ = 
j( 7 r/TO)sgn(x — x'), while all other commutators vanish. 
The charge density is given by p{x) = dxQxjr^ and the 
spin density by s^ix) = dxfiij'x. Here, we measure 
charge (spin) in units of the quasi-particle charge e/m 
(of the electron spin h/2). The non-quadratic parts of 
the Hamiltonian become 

Hx = 4Aa; cos(0i) cos(m0j) (9) 

Hmod = 6:p,cos{9i) cos{m(j)i — mfi). ( 10 ) 

Again, we will focus on the interface between the CDW 
dominated region A^: 7 ^ 0 for x > 0 and the mag¬ 
netic field dominated region p 0 for x < 0. Non¬ 
quadratic terms relevant in the renormalization group 
sense [75] lead to the chiral field 9i being gapped uni¬ 
formly throughout the system, say, 9i = ttM. Pinning of 
other fields is chosen in such a way that the total energy 
is minimized, so 

9i =—{M+ l + 2h), x>0, ( 11 ) 

m 

(j)\ = (j) -\ -(M “h 1 + 2/), X < 0, (12) 
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where M, h, and I are integer-valued operators. The only 
non-trivial commutation relation between them is [h, /] = 
imjArK. The corresponding zero-energy parafermion op¬ 
erator isnnHis] is given by 

a = e*^("+'\ a™ = l. (13) 

The so found ground state is m-fold degenerate and the 
bound states are Z^-parafermions obeying non-Abelian 
braiding statistics [37H1S]. This m-fold degeneracy can 
be explained following Refs. [38H4T] . Let us assume 
that we have a second interface at x = L such that 
= 0 and /2 7^ 0 for a; > L. The spin located in the 
gapped CDW dominated region (0 < cc < L) is given by 
(sz) = 2(Z< — ly)/m, assuming m distinct values, where 

is the quantum number I right/left to this region. By 
analogy, the magnetic field dominated region is charac¬ 
terized by the difference in the charge density between 
the two edges. The parafermion relation is written as 
a<a> = 

Similarly, the interface between the tunneling 
and field {Hx) dominated regions can also host fractional 
charges considered above. We note that Htun rewritten 
in terms of chiral fields becomes of the same form as 
Eq. (|^ but with 0 = 0. 

Exotic superconductivity. Next we consider a TIC in 
proximity to a bulk s-wave superconductor with phase 
0SC and work in the spin filtering regime, see Fig. The 
only proximity-induced superconducting term that can 
open a gap in the spectrum is of the form 

i/L = + H.C., (14) 

where A(,j, is the proximity gap in the TIC. This term 
describes the coherent tunneling of two Cooper pairs out 
of the condensate of the bulk superconductor into the 
TIC edge states. The bosonized version of this term 
reads = 2A(,^cos(20i -|- (j)sc) and has a minimum 
at 01 = 7r(n -|- 1/2) — 0sc for A),^, > 0 and 0i = 7rn — 0sc 
for A),^ < 0, where n is an integer. If it is a relevant term 
in the renormalization group sense m, it opens a par¬ 
tial gap in the spectrum, i.e., the charge degrees of free¬ 
dom are gapped out but the spin degrees of freedom stay 
gapless. The interface between two such regions forms 
a Josephson junction with a Josephson current of Sirm- 
periodicity as a consequence of the dm-fold degeneracy 
of the ground state. 

Opening gaps by interactions. Finally, we comment on 
the possibility of opening gaps in the spectrum solely via 
interactions. For example, the back-scattering exchange 
term Hint — Ai„t(i?|Li)(i?jLi) -|- H.c., which reads in 
bosonized form, 

H^nt = 2Ai„4 cos(26»i), (15) 

opens a gap in the charge sector. This exchange term 
requires overlap between the TIC edge states. The inter¬ 
face between the region with Aint > 0 [0i = 7r(n -|-1/2)] 


and the one with A^t < 0 [di = 7rn] hosts a quasiparticle 
with fractional charge e/2TO. The domain wall between 
and Hint hosts bound states in the charge sector. By 
analogy with the degenerate states considered above [see 
Eqs. ( [TT| ) - (|l^], the ground state in the charge sector is 
2m-fold degenerate. 

The time-reversal invariant two-particle back- 
scattering term described in Refs. m and m can 
result in the opening of a gap of the TI edges and is 
given by H^m = J dx e“^*^^'^0|(i90|)(5'0i)0i -I- H.c. In 
the bosonized form the term becomes 


Hum = 4:Aum cos(26>i) cos(2m6»i). 


(16) 


In contrast to Hint, this term opens a full gap in the spec¬ 
trum with the quasiparticle charge in the system given 
by e/2m with 0i = -kM /2 and 9i = tt{M -I- 1 -I- 2n)l2m. 
If the eigenvalue M is an even (odd) number, {sy) ^ 0 
{{sy) = 0) with {sx) = 0 for all values of M. The in¬ 
terface between Hum and Hmod also hosts bound states. 
The charge field is pinned uniformly with 9i = ttM, 0j is 


pinned according to Eq. (121 and = 7r(h -|- 1 /2 )/to re¬ 
places Eq. Again, two spin fields on different sides 

of the domain wall do not commute with each other, re¬ 
sulting in a degenerate ground state with non-Abelian 
statistics described by parafermion operators. 

Conclusions. We have considered constrictions in two- 
dimensional topological insulators. First, we show that 
such TICs could be used as spin filters operated solely by 
electric fields, i.e. by tuning the chemical potential. Sec¬ 
ond, we demonstrated that the proposed setup could be 
used to generate degenerate fractional bound states with 
non-Abelian statistics. The domain walls occur at the 
interfaces between regions of different gap-opening mech¬ 
anisms. For example, gaps can be opened by magnetic 
fields, tunneling between edges, charge density waves, or 
solely by interactions. We finally note that the proposed 
coupled edge states could be realized not only in TIs but 
also in systems of coupled wires [S3] in the framework of 
strip of stripes [531163] . which is especially relevant for 
the fractional TI regime. 

This work is supported by the Swiss NSF and NCCR 
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